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Abstract — We extend the stabilizer formalism to a class of non- 
additive quantum codes which are constructed from non-linear 
classical codes. As an example, we present infinite families of non- 
additive codes which are derived from Goethals and Preparata 
codes. 

I. Introduction 

Recently, several new non-additive quantum error-correcting 
codes (QECCs) have been constructed that have higher dimen- 
sion than additive QECCs with the same length and minimum 
distance [5], [16], [17]. The first example of such a code is the 
code ((5,6,2)) of [11] which has been found via numerical 
optimization. Afterwards, the code has been identified as the 
span of a particular state and its image under five unitary 
transformations (see also [7]). The recently discovered codes 
((9,12,3)) and ((10,24,3)) (see [16], [17]) start with a so- 
called graph state which corresponds to a stabilizer state, i. e., 
a stabilizer code with parameters [[n, 0, d]] (see [8], [13]). A 
basis of the quantum code is obtained by this initial state 
together with its image under some tensor products of Pauli 
matrices and identity (Pauli operators). The distance between 
any pair of these states can be defined as the minimal weight 
of a Pauli operator transforming one state into the other (see 
below). The problem of finding a code of high dimension 
can be stated as finding a maximal clique in a search graph 
whose vertices are all images of the initial state. There is 
an edge between two states if their distance is at least the 
prescribed minimum distance. Using the formalism of graph 
states. Cross et al. show in [5] that it is sufficient to consider 
Z-only operators in order to define the search graph, but this 
has the disadvantage that the distance between two of these 
states is not necessarily equal to the number of Z operators in 
the tensor product. Moreover, constructing a code {{n,K,d)) 
requires to find a clique of size K in a search graph with 
2" vertices. Fixing one basis state, the graph can be slightly 
simplified. However, for the code ((10,24,3)) the simplified 
graph still has 678 vertices and 149.178 edges. 

A different approach for constructing non-additive QECCs 
based on Boolean functions and projection operators has 
been presented in [1]. Evaluating the Boolean function at the 
projection operators, the code is given as the sum of the image 
of products of the projections. Finally, we mention the non- 
additive QECCs obtained by the method given in [12]. Those 



codes have the same parameters as the additive CSS codes 
that can be obtained using the same underlying binary codes. 

In this paper, we extend the approach of combining unitary 
images of stabilizer states to arbitrary stabilizer codes as 
starting point. For a stabilizer code [[n, k, d]], the search graph 
has only 2"^*^ vertices, and a clique of size K yields a 
quantum code of dimension A' x 2*^. What is more, we present 
some infinite families of non-additive quantum codes which 
are constructed using non-linear binary codes, avoiding the 
NP-hard problem of finding a maximal clique in the search 
graph. Finally, we show how to obtain encoding circuits for 
the resulting non-additive codes using encoding circuits for 
the underlying stabilizer codes. 

II. A BRIEF REVIEW OF THE STABILIZER FORMALISM 

We start with a brief review of the stabilizer formalism for 
quantum error-correcting codes and the connection to additive 
codes over GF{4) (see, e.g., [4], [6]). A stabilizer code 
encoding k qubits into n qubits having minimum distance d, 
denoted by C = [[n, /c,(i]], is a subspace of dimension 2'^ of 
the complex Hilbert space (C^)®" of dimension 2". The code 
is the joint eigenspace of a set of n — A: commuting operators 
Si, . . . , Sn-k which are tensor products of the Pauli matrices 

^-^(i o)' ^^^(* o')' ^^^(o -l)' 

or identity. The operators Si generate an Abelian group S 
with 2"^'"' elements, called the stabilizer of the code. It is 
a subgroup of the n-qubit Pauli group Vn which itself is 
generated by the tensor product of n Pauli matrices and 
identity. We further require that S does not contain any non- 
trivial multiple of identity. The normalizer of S in Vn, denoted 
by J\f, acts on the code C ~ [[n, k, d]]. It is possible to identify 
2k logical operators Xi, . . . , and Zi, . . . ,Zk such that 
these operators commute with any element in the stabilizer S, 
and such that together with S they generate the normalizer 
J\f of the code. The operators Xi mutually commute, and so 
do the operators Zj. The operator Xi anti-commutes with the 
operator Zj if i = j and otherwise commutes with it. 

It has been shown that the n-qubit Pauli group corresponds 
to a symplectic geometry and that one can reduce the problem 
of constructing stabilizer codes to finding additive codes over 



GF{4) that are self-orthogonal with respect to a symplectic 
inner product [3], [4]. Up to a scalar multiple, the elements 
of Vi can be expressed as a^a^ where (a, b) e F2 is a 
binary vector. Choosing the basis of GF{A), where 

w is a primitive element of GF{4) with lu^ + oj + 1 ~ 0, we 
get the following correspondence between the Pauli matrices, 
elements of GF(4), and binary vectors of length two: 
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This mapping extends naturally to tensor products of n Pauli 
matrices being mapped to vectors of length n over GF{4) or 
binary vectors of length 2n. We rearrange the latter in such a 
way that the first n coordinates correspond to the exponents 
of the operators cr^. and write the vector as (a|6), i.e., 

g = a^.^al^ ® . . . ® a^a^ = (a|fe) = (.g^'l/). (1) 

Two operators corresponding to the binary vectors {a\b) and 
(c|c?) commute if and only if the symplectic inner product 
a ■ d ~ b ■ c = 0. In terms of the binary representation, 
the stabilizer corresponds to a binary code G which is self- 
orthogonal with respect to this symplectic inner product, and 
the normalizer corresponds to the symplectic dual code G* . 
The stabilizer together with the logical operators Zj generate a 
self-dual code. In terms of the correspondence to vectors over 
GF(4), the stabilizer and normalizer correspond to an additive 
code over Gi^(4) and its dual with respect to an symplectic 
inner product, respectively, which we will also denote by C 
and C* . The minimum distance d of the quantum code is 
given as the minimum weight in the set G* \ C which is 
lower bounded by the minimum distance d* of G* . If d — d* , 
the code is said to be pure, and for d> d* , the code is said 
to be pure up to d* . 

III. The union of stabilizer codes 

Note that we have defined a stabilizer code C as the joint 
eigenspace of the commuting operators Si generating the 
stabilizer S. The term stabilizer suggests that the code is 
the joint +1 eigenspace of the operators. However, for each 
of the generators Si we may choose either the eigenspace 
with eigenvalue +1 or the eigenspace with eigenvalue — 1. 
This gives rise to a decomposition of the space (C^)*^" into 
2n-fe mutually orthogonal spaces which can be labeled by the 
eigenvalues of the n~ k generators Si, or equivalently, by the 
characters x of the stabilizer group S. Moreover, the n-qubit 
Pauli group Vn acts transitively on these spaces. 

From now on we fix the code C as the the joint +1 
eigenspace corresponding to the trivial character. Let t E Vn 
be an arbitrary n-qubit PauU operator. Then we can define a 
character of S on the generators Si as follows 

I +1 if t and Si commute, 
Xt{Sij := < 

1—1 if t and Si anti-commute. 



As the elements of the normalizer J\f commute with all 
elements of the stabilizer S, two elements <i and t2 define 
the same character if t2^ti £ Af. Hence the set of characters 
corresponds to the cosets of Af in Vn- If T is a set of coset 
representatives, we can write the decomposition of the full 
space as 

(C2)«5»^0tC. (2) 

teT 

Note that measuring the eigenvalues of the generators 5*^ 
projects onto one of these space tC, corresponding to the 
character xt given by the sequence of eigenvalues. In terms 
of the classical codes, the eigenvalues correspond to an error- 
syndrome which is obtained by computing the symplectic 
inner product of the received vector with the n — k vectors 
corresponding to the generators of the stabilizer, i. e., a basis of 
the code G. For all vectors of the dual code G* corresponding 
to the normalizer J\f, the inner product is zero. So the different 
spaces tC correspond to cosets G* + 1 of the code G* . 

As for a fixed code C two spaces tiC and are either 
identical or orthogonal, we can define the distance of them as 
follows: 

dist(tiC, taC) min{wgt(p) : p e V,, \ phC = tjC}. (3) 

Here wgt(p) is the number of tensor factors in the n-qubit 
Pauli operator p that are different from identity. Clearly, 
dist{tiC ,t2C) = dist(i2"^tiC,C). The distance ^ can also 
be expressed in terms of the associated vectors over GF(4). 

Lemma 1: The distance of the spaces tiC and equals 
the minimum weight in the coset G* + ii — 12, where we use ti 
to denote both an n-qubit Pauli operator and the corresponding 
vector over Gi^(4). 

Proof: Direct computation shows 

dist(iiC, t2C) = dist(G* + ti, G* + t2) 

= dist(G* + (ti-t2),G*) 

= min{wgt(c + ii -12) : c e G*} 

= min{wgt(i;) : v £ G* + ti — 12}- g 

With this preparation, we are ready to present the general 
construction of the union of stabilizer codes (see also [7]). 
The quantum code will be defined as the span of some of the 
summands in (|2]l. 

Definition 2 (union stabilizer code): Let Cq = [[n, fc.c?o]] 
be a stabilizer code and let Tq = {ti, . . . , Ik} be a subset 
of the coset representatives of the normalizer Ao of the code 
Co in Vn- Then the union stabilizer code is defined as 

C = iCo. 

teTa 

Without loss of generality we assume that To contains identity. 
With the union stabilizer code C we associate the (in general 
non-additive) union normalizer code given by 

G* ^ U Go + i = {c + : c e Go*, i = 1, . . . , A'}, 
teTo 

where Gq denotes the additive code associated with the 
normalizer Aq of the stabilizer code Co. We will refer to 





5f 


\ 




aZ 

^n-k 










Z^k 


zl 












xl 


J 












j 



Fig. 1. Arrangements of the vectors associated witli a union stabilizer code. 

both, the vectors ti and the corresponding unitary operators, 
as translations. 

A union stabiHzer code can be defined in terms of binary 
vectors as shown in Fig.fT] The first n — k rows correspond to 
the binary vectors (cf. ([TO associated with the generators 5^ 
of the stabilizer S of the code Cq. They generate the classical 
code Cq. The next k rows correspond to the logical operators 
Zj, followed by the k logical operators Xi. The last K rows 
correspond to the K translations defining the cosets of the 
classical code Cq and the unitary images of the stabilizer 
code Co, respectively. We use curly brackets to stress the fact 
that the set of operators Tq need not be closed under group 
operation. In general, the quantum code is not invariant under 
these generalized logical X-operators. On the other hand, if 
Tf) is closed under group operation, the resulting code will be 
a stabilizer code where a basis of Tq defines an additional set 
of logical X-operators. 

Theorem 3: Let C be a union stabilizer code as in Definition 
[2] The dimension of C is |7o|2'' — K2^, and the minimum 
distance is lower bounded by the minimum distance d of the 
union normalizer code C* . 

Proof: As 7^^ is a subset of the coset representatives 
of the normalizer Ao, the spaces tiC^, each of which has 
dimension 2^, are mutually orthogonal. Hence the dimension 
of the union code is K2^ . Fixing an orthonormal basis {\cj) : 
j = 1, . . . , 2'^} of the stabilizer code Co, the set {ti\cj) : i — 
1, . . . , K, j = 1, . . . , 2^} is an orthonormal basis of the union 
stabilizer code. Let E E Vn be an n-qubit Pauli error of weight 
< wgt(£^) < d. For basis states \cij) — ti\cj) E tiCa and 
\ci\j') — ti'\cj') E tiiCo we consider the inner product 

{c^.J\E\ce^r)^(c,\tlEU,\c,,). (4) 

For i ^ i', we have dist(tiCo, ti'Co) — min{wgt(c + ti — t^') : 
c G Cq} > d, and hence Q vanishes for wgt{E) < d. For 



i — i', we get 

{c,\tlEU\c,,)^±(cj\E\c,,). (5) 

As the code Cq is pure up to the minimum distance of Cq C 
C*, equation (j5]l vanishes as well. ■ 
Remark 4: We note that similar to stabilizer codes, the true 
minimum distance of a union stabilizer code might be higher. 
The true minimum distance is given by 

min{dist(c + ti,c + U') : ti, ti' E %, 

c,c' EC*\c + t,-ic' + t,,)(^Co} 
= min{wgt(w) : V E (C* - C*) \ Cq}, 

where (C* — C*) = {a — b : a,b E^C*} denotes the set of all 
differences of vectors in Cq, and Co is the symplectic dual of 
the additive closure of C*. 

In order to construct union quantum codes, we may start with 
a stabilizer code Cq and use a search graph whose vertices 
are the mutually orthogonal translates {tCo : t E T^} of the 
stabilizer code. Two vertices are connected by an edge if and 
only if the distance between them is at least d, where d is 
the desired minimum distance. For simplicity, we also require 
that the code Co is pure up to do > d. The distance between 
two translates can be computed using Lemma [T] This allows 
to use stabilizer codes Co of arbitrary dimension, and hence 
allows to go beyond the case of stabilizer states (or graph 
states) as, e.g., in [5], [17]. We note that the construction of 
non-additive quantum codes of [2] is also based on taking the 
union of orthogonal images of a stabilizer code. 

IV. Union Stabilizer Codes from Binary Codes 

A. CSS-like codes 

Given two linear binary codes Ci = [n,ki,di] and C2 = 
[n, ^2,^2] with C^ C Ci, the so-called CSS construction 
(see, e.g., [14]) yields a quantum error-correcting code C — 
[[n, fci + fc2 — n, d]] with d > min(di, ^2)- Starting with this 
CSS code, we consider unions of cosets of the binary codes 
Ci, i. e., 

a= \J c, + t« 

such that the minimum distance of the codes Ci is at least 
d < d. Using the translations {(t^^^l^^^)) ; ^(i) ^ Ti,t^^^ E 
7^} we obtain a CSS-like union stabilizer code of dimension 
\Ti\ ■ |72| • 2'^i+'^2^" whose minimum distance is at least d. 
If Ci = iuii) ^2 = ((^) generator matrices of 
the codes Ci, where Hi is a generator matrix of the dual code 
Cj^, the corresponding vectors are as shown in Fig. [2] 

B. Enlargement construction 

Steane has presented a construction that allows to increase 
the dimension of a CSS code [14]. For this, he starts with the 
CSS construction applied to a binary code C = [n, k, d] which 
contains its dual, yielding a CSS code Cq = [[n, 2A: — n,d]]. 
Using a code C = [n,k' > k + l,d'] which contains C, he 
obtains a quantum code [[n, k + k' — n,mm{d, \3d' /2'])]]. The 
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Fig. 2. AiTangements of the vectors associated witli a CSS-like union 
stabilizer code. 



resulting code can also be considered as a union stabilizer 
code. If is a generator matrix of the complement of C 
in C and A is a fixed point free, invertible linear map, the 
translations can be defined as 

To = {{vD\vAD) : v e F^'-*^}. 

The key observation [14] for proving the lower bound on 
the minimum distance is that the weight of an operator g = 
{g'^\g^) can be expressed in terms of the Hamming weight 
of the binary vectors and their sum: 

wgt((5^|5^)) = ^(wgt(5^) + wgt(g^) + wgt(g^ + g^)). 

As To is closed under addition and the properties of A ensure 
that Q^t^^t^^Q for any non-zero element {t^\t^) G %, 
the weight of all three binary vectors is lower bounded by d! . 

V. Quantum Codes from Reed-Muller, Goethals, 

AND PREPARATA CODES 

Using the CSS-like construction of the previous section, we 
now construct some families of non-additive quantum codes. 
For this, we use the Goethals codes Q{m) and the Preparata 
codes which are nonlinear binary codes of length n ~ 2™ for 
m > 4, m even. Some of the properties of these codes are 
summarized as follows (see, e. g., [10]): 

• Both the Goethals code Q{m) and the Preparata code 
V{m) are unions of cosets of the Reed-Muller code 
Tl{m) := RM(m — 3, m). Furthermore, they are nested 
subcodes of RM(to — 2, m), i. e., 

RM(m - 3, m) C g{m) C 'P(m) C RM(to - 2, m). 

• The parameters of the codes are 

RM(m - 3, m) = n{m) = [2™, 2" - (™) - m - 1, 8] 
g{m) = (2",22'"-3m+i^8) 

V{m) = (2",22'"-2m^6) 
RM(m-2,m) = [2", 2" - m - 1, 4]. 
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Fig. 3. Transformation of the union stabilizer code given by the inverse 
encoding circuits Qi and Qc- 



Steane has constructed a family of additive quantum codes 
from Reed-Muller codes [15]. The codes are obtained applying 
the enlargement construction of [14] to the chain of codes 

RM(r, m) C RM(r, m)^ = RM{m - r - 1, to) 
C RM(to - r,m). 

In particular, for r = 2 and to > 5 this yields additive QECCs 
Ci = [[2™, 2™ - (™) - 2to - 2, 6]], while using only the CSS 
construction, one obtains Co = [[2™, 2™ - 2(™) - 2to - 2, 8]]. 

As the Goethals code Q (to) is the union 
cosets of Tl{m), we can construct a CSS-like union stabilizer 
code based on Cq. The minimum distance of the resulting 
non-additive code is 8 and its dimension is KQdhTi{Co) = 

Replacing the Goethals code by the Preparata code V{m), 
we have Kp = 2(2)^'"+^ cosets of TZ{m). This results in a 
CSS-like union stabilizer code with minimum distance 6 and 
dimension ii:|>dim(Co) = 2"^"^-^""'. 

Both the union stabilizer codes based on Goethals codes 
and those based on Preparata codes are superior to the additive 
codes derived from Reed-Muller codes. The parameters of the 
first codes in these families are as follows: 
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However, applying the enlargement construction to extended 
primitive BCH codes results in stabilizer codes with parame- 
ters [[2™, 2" - 5to - 2, 8]] and [[2™, 2" - 3to - 2, 6]] (see 
[14]). 

VI. Encoding Circuits 

In [9] we have shown how to compute a quantum circuit 
consisting of Clifford gates only that transforms any stabilizer 
S given by the binary {n ~ k) x 2n matrix {X\Z) into the 
stabilizer of a trivial code given by (0|/0), where / is an 
identity matrix of size n — k. The corresponding trivial code 
corresponds to the mapping [(/)) i~> |0 ... 0) [(/)). We denote 
the resulting quantum circuit that corresponds to the inverse 
encoding circuit of {X\Z) by Qi. Note that we can apply Qi 
to all the operators defining the code as illustrated in Fig. |3] 
Further note, that for the trivial stabilizer code, the "encoded" 
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Fig. 4. Inverse encoding circuit for the non-additive code ((5, 6, 2)). Tlie first set of gates including the Hadamard transformations implements the inverse 
encoding circuit Q\ for the stabilizer code [[5, 0, 3]], followed by 5 CNOT and 2 Toffoli gates implementing the classical circuit Qc- 



X- and Z-operators are weight-one Pauli operators cTj. and a^, 
respectively, acting on the last k qubits. As the transformed 
translations define cosets of the normalizer, we may 

choose them such that they are tensor products of operators 
and identity acting on the first n — k qubits only. Then we have 
the trivial union code spanned by a set of K2^ basis states 
of the form |ci)|j), where : j = 0, . . . , 2*^' — 1} is the 
computational basis of k qubits and {ci : i = 0, . . . , K — 1} 
is a set of bit strings of length n — fc. In order to obtain a 
standard basis for our input space of dimension K2'', we need 
a quantum circuit Qc mapping \ci) i-^ \i) for i = 0, . . . , A' — 1. 
Note that this is a purely classical circuit which can be 
realized, e. g., using ax, CNOT gates, and Toffoli gates. 

We illustrate this for the non-additive code ((5, 6, 2)) which 
is a union stabilizer code derived from the stabilizer state 
Co = [[5,0,3]]. For a stabilizer code with k = 0, there are 
no encoded X- and Z-operators. So the code ((5,6,2)) is 
specified by five generators of the stabilizer and six transla- 
tions. Using an inverse encoding circuit Qi, these operators 
are transformed as follows: 
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It remains to find a (classical) quantum circuit Qc that maps 
the six binary strings on the right hand side of (|6]l to say 
the binary representations of i = 0, . . . , 5. Using a breath-first 
search among all circuits composed of cr^, CNOT, and Toffoli 
gates, we found the minimal realization shown together with 
the circuit Qi in Fig. |4] 

VII. Conclusions 

The approach presented in this paper generalizes naturally 
to the construction of non-additive quantum codes for higher 
dimensional systems. In order to obtain other families of non- 
additive quantum codes, it is interesting to study classical non- 
linear codes which can be decomposed into cosets of linear 
codes, similar to the Preparata and Goethals codes. 
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